Decoherence and Initial Correlations in Quantum Brownian Motion 
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We analyze the evolution of a quantum Brownian particle starting from an initial state 
that contains correlations between this system and its environment. Using a path integral 
approach, we obtain a master equation for the reduced density matrix of the system finding 
relatively simple expressions for its time dependent coefficients. We examine the evolution of 
delocalized initial states (Schrodinger's cats) investigating the effectiveness of the decoherence 
process. Analytic results are obtained for an ohmic environment (Drude's model) at zero 
temperature. 
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I. INTRODUCTION 

The study of the decoherence process has received in- 
creasing attention in recent years Q. In fact, it has been 
recognized that decoherence is of fundamental impor- 
tance in understanding the nature of the fuzzy boundary 
between the quantum and the classical domains. The 
nature of this boundary has been under scrutiny both 
from the theoretical and from the experimental point of 
view m . The basic physics of decoherence is very sim- 
ple: interaction with the environment tends to prevent 
the stable existence of the vast majority of the states 
of the Hilbert space of macroscopic quantum systems. 
Thus, coherent superpositions of macroscopically distinct 
states tend to decay very rapidly (on a short decoherent 
time-scale) into mixtures preventing the observation of 
delocalized (Schrodinger's cat) states. In some sense, the 
interaction with the environment enforces the existence 
of "environment induced" super-selection rules selecting 
the very few states in which classical systems are found. 
To answer specific questions concerning the effectiveness 
of the decoherence process for inducing classical behav- 
ior in a particular system one certainly has to analyze 
detailed models describing the actual physical situation. 
However, some generic features of decoherence have been 
analyzed for classes of models which appear in a variety 
of physical circumstances. The paradigmatic model for 
such studies, which we will reanalyze in this paper, has 
been the linear quantum Brownian motion (QBM) which 
is characterized in the following way: A Brownian parti- 
cle (whose coordinate we denote with q) evolves in one 
dimension while interacting with an environment formed 
by a collection of independent harmonic oscillators (with 



coordinates £„). The Lagrangian of the system environ- 
ment ensemble is: 



L s (q) 



L s (q) + L se (q,£) 
\k 2 V(q) 
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(3) 



Notice that, for convenience, L se contains both the free 
Lagrangian for the environment as well as the interaction 
term (including a potential frequency renormalization). 
For the potential we will simply consider V(q) — uj 2 ,q 2 /2. 

This model is such that only a few parameters are im- 
portant to describe the effects the environment produce 
on the system. One of such "parameters" is the so called 
spectral density of the environment which measures the 
density of oscillators with any given frequency and the 
strength with which they couple to the system. This 
function is defined as 



6{lo - Lo n ). 



(4) 
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The other important ingredient required to determine the 
effect of the environment is the initial state. The simplest 
such state for the system-environment ensemble is a fac- 
torizable state where the total density matrix is just a 
product: p se — p s ® p e - For this type of states, and un- 
der a variety of assumptions, the decoherence process has 
been analyzed for QBM models . In this 

paper we will study the evolution of the reduced den- 
sity matrix and, in particular, analyze the effectiveness 
of the decoherence process for a much wider class of ini- 
tial states than the ones analyzed so far in the literature. 
We will allow for initial states for which the initial density 
matrix is not factorizable (i.e., states containing system- 
environment correlations) . Our work will be based on the 
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use of techniques and results that have been elaborated 
and clearly exposed by Grabert et al in [Q. 

One of the most practical tools for analyzing the evo- 
lution of a quantum open system [jl2 13 is the evolution 
equation for the reduced density matrix. This is known 
as the master equation and its properties for the QBM 
model have been extensively analyzed in the literature 
ji^,[l5]|||| . However, only relatively recently it has been 
realized that the structure of the master equation for lin- 
ear QBM models is always remarkably simple || (see also 
[^,|lTJ). Thus, for the case of factorizable initial states it 
has been shown that the exact master equation for lin- 
ear QBM is always local in time having time dependent 
coefficients. A variety of derivations of this exact mas- 
ter equation, valid for environments with general spectral 
densities in initial states of arbitrary temperature, have 
been given so far in the literature §,|,|Jl7|]l|] . This type 
of equations have been used to analyze a rather wide va- 
riety of problems (see [^9| , p0|j2l| l for references of use of 
master equation and related techniques in context rang- 
ing from cosmology to quantum optics). 

In this paper we will generalize previous work on mas- 
ter equation for QBM models allowing for a more gen- 
eral class of initial states and finding the general form of 
the exact master equation, Our equation reduces to the 
previously known one |6|,|8| for the case of vanishing ini- 
tial correlations. As an aside, we present a very simple 
derivation of the master equation and find rather conve- 
nient and manageable formulae for the time dependent 
coefficients. Using them, we analytically solve a simple, 
but physically relevant, example where all the coefficients 
can be computed (Drude's model of an ohmic environ- 
ment at zero temperature). Of course, this is not the first 
time the QBM model with non-factorizable initial states 
has been analyzed. As we mentioned above, the method 
we apply here has been developed and used by others (see 
JDlOl)- However, to our knowledge, neither the struc- 
ture of the master equation has been investigated in this 
case before, nor the effectiveness of decoherence has been 
examined (except for the work in |22]| which we general- 
ize here). In some sense, our paper is part of an effort 
to relax the usual assumptions behind simple models of 
decoherence (further work towards a more complete "de- 
construction" of decoherence is presented elsewhere [E3| ) . 

The paper is organized as follows: In Section 2 we de- 
scribe the class of initial states we analyze and introduce 
the concept of preparation function. In Section 3 we de- 
scribe the formalism following the scheme presented in 
P|. In Section 4 we obtain the master equation de- 
scribing its properties and studying the behavior of its 
coefficients. In particular, we find analytic expressions 
for Drude's model at zero temperature. In Section 5 
we study the evolution of two types of dclocalized ini- 
tial states (Schrodinger's cat states) consisting of a su- 
perpositions of Gaussian wave-packets. In both cases we 
analyze the evolution of the Wigner function analyzing 



the efficiency of the decoherence process. Finally, in Sec- 
tion 6 we summarize our conclusions. Appendices 1, 2 
and 3 contain useful formulae which we do not include 
in the main text to prevent overloading it with too many 
equations. 



II. INITIAL CONDITIONS 

We are interested in studying the following type of ini- 
tial conditions: the system and the environment have in- 
teracted for a very long time so that they reached an equi- 
librium state represented by the density matrix pp. At 
the initial time (t — 0~) we make a measurement on the 
system only. As every result of the measurement is asso- 
ciated with a projection operator P acting on the Hilbcrt 
space of the system, the state of the system+environment 
ensemble after such ideal measurement is 



Po 



Tr(Ppp) 



(5) 



It is clear that the above is not a product state since 
it contains correlations between the system and the en- 
vironment, which are inherited from the ones already 
present in the pre-measurement thermal state pp. There- 
fore the usual techniques are not applicable for describing 
the evolution of state (JsJ) . In what follows we will present 
a method enabling us to study the fate of a whole class of 
states which includes (||) as a particular case. In general, 
the initial states we will consider are of the form: 



Po = A jPP A 'j 



(6) 



where A, and A'j are operators (not necessarily projec- 
tors) acting on the Hilbert space of the system. For 
equation (g) to represent the state following a perfect 
measurement on the system, the above sum must col- 
lapse onto a single term and Aj = Aj = P/(Tr(Ppp)) 1 / 2 
where P is a projector. 

At this point it is convenient to introduce some nota- 
tion: It will turn out to be useful to describe the initial 
state of the system in terms of a "preparation function" 
\{q, q, q' , q'). This function, which parametrizes the devi- 
ation of the initial reduced density matrix of the system 
from its thermal equilibrium form, is defined in terms of 
the matrix elements of the operators Aj,A'j as: 



A(g, q, q', (?) 
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< q\A 3 \q> < q'\A'Aq' > 



(7) 



Using this definition, it is easy to show that the total 
density matrix in the coordinate representation is: 

Po(q,£,q',0 = / dq dq X(q,q,q',q) pp{q, q , £')■ (8) 
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and that the initial reduced density matrix is 

P(<l,<l')= / dq dcf\(q,q,q',q') pp{q,q'). (9) 



III. EVOLUTION OPERATOR 

As described by Grabert and others |ll]], it is possible 
to modify the usual Feynman- Vernon method p4| (which 
is only applicable for factorizable initial states, see gj) 
to study the evolution of initial states of the form (q). 
The main results, whose derivation we will briefly outline 
below are the following (see Jll| for more details): 

(i) The reduced density matrix of the system can be 
obtained by "evolving" the preparation function in the 
following way: 

P(q, q', t) = J dqi dq'i dq dq'J{q, q' , t, q^q-, q, q') 

xA(«,9,i?)> ( 10 ) 

(ii) The evolution operator J has a simple represen- 
tation as a triple path integral over trajectories of the 
system. In this path integral representation the effect 
of the environment is present through a Generalized In- 
fluence Functional which provides a nontrivial weight to 
triplets of system's trajectories (see below), 

(iii) For general linear models (i.e., an environment of 
independent oscillators with an arbitrary spectral den- 
sity) the problem is exactly solvable. Thus, if the ac- 
tion is quadratic in the environmental coordinates, the 
Generalized Influence Functional can be easily computed. 
Moreover, if the theory is linear for the system, a closed 
expression for the evolution operator J can be obtained. 

Let us now describe how to demonstrate results (i)- 
(iii) quoted above. The validity of (i) can be simply seen 
by writing the full density matrix at time t in terms of 
the initial density matrix as 



P(q,£,q',£,',t) = J dKdq'd^p^^q',^) x 

x K(q,lt, qi ,^) K*(q',£',t,ql,&. (11) 

where, K is the evolution operator of the complete wave- 
function. Equation ( [To| ) is obtained by expressing the ini- 
tial density matrix in terms of the preparation function 
(as in (||)) and by tracing over the environment coordi- 
nates £. Doing this we also obtain the explicit form of 
the evolution operator J: 



J(q,q',t,q t ,q' l ,q,q') = J d&d&d^f p p (q, q' , £■) x 

x K{q,^,t,qiAi) K*tf,Zf,t,&&) ■ (12) 

To find a simple path integral representation for this 
evolution operator (property (ii) above) we can first ex- 
press the full evolution operator K as a sum over histories 
of the system and the environment: 



K{q,^t;q h & = J DqDi 



(13) 



where the integration paths must satisfy the boundary 
conditions: 

q(0) = qi, q(t)=q, £(0) = £(t) = £. (14) 

Replacing this into ( |l2|) and expressing the total action 
as a sum of free and interaction terms, the evolution op- 
erator can be written as a path integral: 



J{q,j,t,<k,&4,?) = J D q D q' e^M-^te] 
x J d£ f d£i dZl /£>/?(<?, Cg'.f) 

x / D£D? e iS "[^- <s " I 9 ''*'!, (15) 



This formula is not of the desired form yet: we have a 
double path integral over the system's trajectories with 
an integrand which is not only a functional of these tra- 
jectories but also a function of q and q'. To overcome 
this difficulty (which comes precisely from the fact that 
the initial state contains correlations, which are present 
in the thermal density matrix pp) we can use the Eu- 
clidean path integral representation for a thermal equi- 
librium density matrix. Thus, matrix elements of pp can 
be written as: 



Dq D£ e" 



(16) 



where the integral is over Euclidean paths satisfying the 
boundary conditions: 



2(0) 



Using this, we can write the evolution operator as the 
following triple path integral (property (ii) above): 

J(q, q\ t, q^q^q, q') = J DqDq' Dq ^M-^]-^m 

xF[q,q',q], (18) 

where the "Generalized Influence Functional" F[q,q',q] 
is defined as: 



F[q,q',q] 



iS. e -iS'-S* 



(19) 



The boundary conditions on the path integral over the 
environment histories are such that all these integrals are 
tied together: final conditions for £ and £' coincide (be- 
cause of the final trace over the environment) while their 
initial conditions are connected via the Euclidean trajec- 
tories. For this reason, the above integral is denoted as 
a functional integral over a "closed time path" pa]. 
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The above considerations are applicable for arbitrary 
interactions. From now on we will restrict ourselves to 
discuss linear QBM models which have the virtue of being 
explicitly solvable allowing an explicit y calculation of the 
propagator J. To do this, one can first notice that, if the 
system-environment interaction is linear in the environ- 
mental coordinates £, the generalized influence functional 
can be exactly computed since the path integral ([l9| ) is 
Gaussian. For the bilinear interaction given by p) the 
result for the "Generalized Influence Action" q' , q] is 
(where F[q, q'q] = exp i$[q, q'q]) is: 



<I>fa;, r 



Z S 

q]= l jdsjdu v{ s-u) x(s) x(u) 





t 



ds / du r/(s — u) x(s) r(u) 



o o 
t 



— ri ds r](s) x(s) 



p p 

J dr J da k(r — a) q(r) q'(cr) 




P t 

+ J dr J ds k*{s - ir) q(r) x{s), (20) 
o o 

where, for convenience we used "sum and difference" co- 
ordinates defined as 



x = q-q 



q + q' 



(21) 



The kernels appearing in the Influence Action (g0|) are 
determined by the spectral density and the initial tem- 
perature 1/(3 (see Appendix 1 for the explicit form of 
these kernels). Here, we would just like to mention that 
the first two lines of equation ( pp| ) contain the usual re- 
sult derived in the absence of initial correlations. Thus, 
kernels v(s) and rj(s) are, respectively, the noise and dis- 
sipation kernels: 



/duu 
— I(lo) coth/3o;/2 cosws 
7T 

oo 

j 



2 ±^ 



71 — - — OO ' 



uj 2 + vl 



COS UJS, 



(22) 



dw 



-I(lu)cos(ws), 



between real and Euclidean trajectories in the influence 
functional. As will be seen below, one of the effects of this 
coupling is to produce an "effective noise kernel" R(s, u) 
(which determines the strength of the diffusive effects). 
The effective noise kernel is not homogeneous in time and 
can be written as: 



i?(s, u) = v{s — u) + V C orr{s, u) 



(24) 



where the explicit form of the "correlational noise" v corr 
is given in Appendix 1. 

Computing the propagator J is straightforward if the 
action is quadratic in the system's coordinates. The final 
result of this simple but tedious calculation (see details 
in 11 1) is: The propagator is a Gaussian function, 



J(x,r,t,Xi,n,x,f) = a ^M, r ^,f) ^ 



(25) 



where the exponent S is a second degree polynomial of 
its arguments which reads: 

S(x, r, t , Xi,ri, x,f) — i {ct\f 2 + a^x 2 ) 

+ a-zixiTi + xr) + a4 Xir + 015 xti 

+ «6 Xif + ioi-j XiX + as xr + iag xx 

+ i (a w x 2 + +an SjS + a 12 x 2 ). (26) 

The explicit formulae for the coefficients ao, . . . ,a±2 
are given in Appendix 1. In general, these time depen- 
dent coefficients arc determined by the spectral density 
and the initial temperature. Here, we will just mention 
a few simple properties of the coefficients: 

(i) cko j us t ensures the normalization (preservation of the 
trace of the density matrix) and it is therefore determined 
by the other coefficients (explicitly, a^ = a 2 ai/16ir 3 ). 

(ii) Coefficients ai, ai are time independent. These coef- 
ficients determine the reduced density matrix in thermal 
equilibrium. Thus, if we denote the position and momen- 
tum dispersion (in thermal equilibrium) respectively as 



and 



.c. if 



2 „ 2 ^ 

q/3 =< q >p, 



pI =< p 2 >p, 



(27) 



then a\ = 1/2<7| and c*2 = p^/2. The explicit form of 
q 2 o and p|, which are temperature dependent is given in 
Appendix 1. 

(iii) Some of the coefficients (a^, a± and Q5) only depend 
upon the spectral density of the environment. Explicitly, 
we have: 



G 



1 



«5 



r 2 

4-a>. 



(28) 



(23) where the function G(t) is a solution of: 



where v n are the Matsubara frequencies v n — 2im/ (3. 
The initial correlations are responsible for the coupling 



G(t) + iu 2 G{t) + — j dt'r)(t - t') G{t') = 0, 



(29) 
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satisfying the boundary conditions 

G(0) = 0, G(0) = 1. 



(30) 



(iv) The other coefficients (a&, . . . , au) depend on G(t) 
and on the noise kernels appearing in the influence func- 
tional. Therefore, they are determined by the spectral 
density and the initial temperature of the environment. 
In the absence of initial correlations (i.e., if we disregard 
the interactions between the system and the environment 
in the Euclidean integrals) the coefficients a§, . . . ,a.g are 
identically zero and the propagator does not mix the co- 
ordinates x, f with the rest. 



IV. MASTER EQUATION 

A. A simple derivation 

Knowing the propagator J, it is possible to find the 
evolution equation for the reduced density matrix. This 
is the so-called master equation, which can be easily de- 
rived following a simple method outlined by one of us in 
||. First, we explicitly evaluate the time derivative of 
the evolution operator which, taking into account (25), 
is of the form J = P 2 J where P 2 is a second degree poly- 
nomial in the variables x, r, Xi, f. Next, we multiply 
this expression by the preparation function and integrate 
over the coordinates Xi, x, f. The nontrivial task is to 
rewrite the right hand side of the resulting formula as an 
operator acting on the reduced density matrix. Part of 
this task is simple: terms involving the "final" coordi- 
nates (x, r), can be moved outside the integrals generat- 
ing local terms in the master equation. The real problem 
is to manipulate terms involving the "initial" coordinates 
(which are being integrated out). To do this, as explained 
in ||, we can take advantage of the following identities, 
which can be straightforwardly derived from eqs (Eq) : 



Xi J = — {—id T — a^x) J 

TjJ = - — d x J (a 3 r H — (d r - i a-ix) + 

+ a^f + iagx) J. 



(31) 



(32) 



Using these equations we can eliminate the initial vari- 
ables Xi, Ti in favor of the rest. At this point we may note 
that the existence of initial correlations (reflected in the 
presence of non- vanishing coefficients as, ag) prevents us 
from completely accomplishing our goal since the right 
hand side of ( |32| ) still depends on the integration vari- 
ables f, x. These terms will generate nontrivial contribu- 
tions to the master equation whose form will be described 
below. Using the above formulae, after some algebra we 
obtain the exact master equation for linear QBM, which 
reads 



Kq,j,t) = l Q (d 2 -d 2 ) - \u)*(t){q 2 p( q ,q',t) 

- 7 (t) (q-q') (d q -d q ,)p(q,q',t) 
-£>l(t) (q-q') 2 p(q,q',t) 
+ lD 2 (t) (q-q 1 ) (d q + d q ,)p(q,q',t) 
+ iC X (t) (q-q') pi(q,q',t) 

-iC 2 (t) (q-q') p 2 (q,q',t) (33) 

cnea It is important to stress that the above master equa- 
tion is exact and valid for all spectral densities and initial 
temperatures. The time dependent coefficients appear- 
ing in ( |33"|) are functions of olq, . . . , ot\ 2 . Explicit formulae 
are given below. 

The interpretation of the terms appearing in the mas- 
ter equation is clear. The first line is just Liouville's 
equation with a renormalizcd Hamiltonian. Thus, the 
environment renormalizes the Brownian particle which 
acquires a time dependent frequency. The form of u)(t) 
is relatively simple (see below) and only depends on the 
function G(t). The second line contains a friction term 
with a time dependent damping coefficient 7(f). This co- 
efficient also has a relatively simple form which is again 
determined by the function G(t). The third line corre- 
sponds to a diffusion term and its presence is of impor- 
tance for studies of decoherence. The diffusion coeffi- 
cient depends both on G(t) and on the noise kernel (its 
time dependence will be examined below in an illustra- 
tive example). The fourth line contains an extra diffusion 
term (called anomalous diffusion in B) which has inter- 
esting effects especially at low temperatures. The last 
two lines make the master equation non-homogeneous. 
In fact, these terms are present because of the corre- 
lated nature of the initial state and prevent the r.h.s. of 
the master equation from being written entirely in terms 
of the reduced density matrix. The two density matri- 
ces p\ and pi are obtained by propagating the initial 
states associated with the preparation functions Ai = fX 
and A2 = xX. Taking into account the definition of the 
preparation function, these states have "density matri- 
ces" pi = {q, p} and p 2 — i[q, p\- It is worth noticing that 
the evolution of pi can also be studied with our formalism 
since (apart from not being normalized) they belong to 
the class of initial conditions defined by @. Therefore, 
the evolution equation for pi is also (^3|), with new in- 
homogeneous terms pij. Thus, a hierarchy of equations, 
which are coupled because of the initial correlations, can 
be derived in this way (see |26| for more details). 

The result we just presented has a remarkable property 
that, at first sight, may appear to be rather tantalizing: 
the master equation ( |33"| ) is local in time (disregarding 
the inhomogeneous terms). Locality of the exact master 
equation, which in the absence of initial correlations was 
previously noticed in [p[p^] , is hard to reconcile with the 
intuitive idea one has about the effects a generic envi- 
ronment may produce. In fact, as such generic environ- 
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merit could produce all sort of non-Mar kovian effects, 
one would expect to find non-local integral kernels in 
the master equation (like the ones appearing in equation 
(p9|)). However, as our exact calculation shows, this is 
not the case for the linear QBM model we are consid- 
ering. Thus, linearity imposes an enormous constraint 
forcing the master equation to be local in time. To un- 
derstand this in a simple way we propose the following 
exercise to the reader: Consider the integro-differential 
equation 

t 

/(f) + 0J o 2 f(t) + i / dt'ri(t - t>) f{t') = 0, (34) 
o 

(which is linear but nonlocal in time). We will show that 
this equation can be easily transformed into the following 
local equation with time dependent coefficients: 



/+ 7 (f)/ + W 2 (f)/ = 



(35) 



The coefficients w{t) and 7(f) (which are precisely the 
same ones appearing in the master equation) can be writ- 
ten in terms of a particular solution of (|34| ) (satisfying 
boundary conditions jl^)) as: 



2 1 . GG-G 2 W 
W = GG-G 2 . 



(36) 
(37) 



Notice that W is just the Wronskian between G and G. 

To demonstrate this, we should first notice that the 
space of solutions of equation (^) has dimension two. 
Therefore the general solution can be written as a linear 
combination of two independent solutions. We can use 
G and G (satisfying initial conditions (|l4j) as a conve- 
nient basis and write f(t) = aiG(t) + a 2 f(t) where a%, a 2 
are two independent constants. These constants can be 
expressed in terms of /, /, G and G as: 



Gf-Gf -Gf + Gf 

fl l = TZi a 2 = 



GG-G 2 



GG-G 2 



(38) 



Replacing this in the nonlocal term of equation ( |34| ) 
and using the fact that G and G satisfy the same equa- 
tion, one can show that f Q dt'r\(t — t')f(t!) = 7(f)/(f) + 
(uj 2 (t) -u>l) f(t). This ends our proof. What this shows 
is a rather trivial feature that is frequently forgotten: 
the future evolution of a function /(f) satisfying a lin- 
ear integro-differential equation like ( |3^ ) does not de- 
pends on its entire history. Indeed, its future behavior 
is uniquely determined by the Cauchy data (the value of 
/ and its derivative) together with the time t at which 
these data are given. The only non-Markovian feature 
of the evolution is that it remembers the time! One may 
argue that equation (|3^) is rather useless: Thus, to solve 
it we must first know u>(t) and 7(f), which means that we 



still need to solve the integro-differential equation. This 
is certainly correct. However, ui(t) and j(t) are "univer- 
sal" functions in the sense that they do not depend upon 
the boundary conditions of the problem at hand. Thus, 
while equation ([35|) is equivalent to ( |34| ) it makes evident 
the fact that the behavior of the solutions is not history 
dependent. 

This simple exercise not only shows how to "local- 
ize" equation ( |34| ) but also makes clear why is the mas- 
ter equation for linear QBM local in time. The non- 
Markovian features we expect to see are rather restricted 
by linearity. In the absence of initial correlations, the 
time dependent coefficients play therefore the very im- 
portant role of providing all the memory effects in the 
evolution of the density matrix. The inhomogeneous 
terms are responsible for carrying the effect of initial cor- 
relations on the evolution of the system. 

Finally, we can write down the equation for the Wigner 
function, which is defined in terms of the density matrix 
as: 



W(r,p,t) = J ^-t 



c p(x,r,t). 



(39) 



Applying the master equation (|33|), it is easy to obtain 
the evolution equation for W{r,p) which results: 

W = {H R , W} PB + 7 (f) d p {pW) + D x (f) d 2 W 

- D 2 (t) d p d r W - d p Woi + C 2 (t) d p W o2 (40) 



where the renormalized (time dependent) Hamiltonian is 
Hr = \p 2 + T)Lo 2 {t)q 2 and {,}pb denotes the standard 
Poisson bracket. This equation, which carries the same 
information than the master equation, is useful for ana- 
lyzing some properties of the solution. In particular, as 
we will show below, it makes transparent the role of the 
anomalous diffusion term. 



B. The time dependent coefficients. 

The expressions defining the diffusive coefficients of 
the master equation are rather complicated at first sight. 
However, we have been able to find the following simple 
formulae: 



(41) 
(42) 



D 1 {t)=[d\,+ 1 {t)d t ,+w 2 {t)] U(t,t') 
D 2 (t)= [d 2 t ,+~,(t)d t ,+u\t)] U(t,t') 

where U(t,t') is an auxiliary function defined as 

t t' 

U{t,t') = J ds J du G{t- s)R(s,u)G{t' - u). (43) 
o o 



G 



and U(t,t') denotes the partial derivative with respect 
to t. Remember that the "effective noise kernel" R(s, u) 
is defined as R(s,u) — v(s — u) + v corr (s,u) where in 
v CO rr(s, u) we incorporate contributions to the noise aris- 
ing from the initial correlations (which are typically rele- 
vant only on very short time-scales). The above expres- 
sion is remarkably simple for the case of uncorrelated ini- 
tial conditions (where R is identical to the standard noise 
kernel v{s)). It can be shown using our equations that the 
known expressions for the diffusion coefficients (found, 
for example, in H) are recovered in this limit. However, 



equations ( [fT|]42| ) are substantially simpler than the usual 
ones (and are valid for a wider class of initial conditions). 

The coefficients appearing in the inhomogeneous part 
of the master equation are also rather complicated but 
we were also able to write them in a simple form as: 



Ci (t) = [d 2 t + i{t)d t + uj 2 (t)] Gae 



(44) 



Before analyzing an explicit example, we should point 
out that the diffusive coefficients are entirely determined 
by the position autocorrelation function (in thermal equi- 
librium), 



S(t) = \ <{?(*), g(0)} > p 



(45) 



In fact, using results found in Appendix 2, we can show 
that: 



U(t,t% =t , 
d t U(t,t')\ t=t , 

d t d t ,u(t,t')\ t=t , 
d t d 2 ,u(t,t% =t , 

d 2 U(t,t')\ t=t , 
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(46) 
(47) 
(48) 
(49) 
(50) 



The equilibrium dispersions p| and q 2 are also deter- 
mined by S(t) through the relations: 5(0) = gj and 
5(0) = — Pg. Thus, knowing the position autocorrelation 
function one could, in principle, compute the diffusive 
coefficients. It can also be seen that the inhomogeneous 
terms in the master equation are also determined by S(t) 
(for example, we have a§ — -jr^ ~ 

It is useful to examine the behavior of all the coef- 
ficients in a particular example. We will consider the 
Drude model, which is characterized by the following 
ohmic spectral density: 



Here, w c is a high energy cutoff frequency below which 
the spectral density is approximately linear in w. The 
integro-differential equation ( |29[ ) can be exactly solved 
by using standard Laplace transform techniques. In this 
way, we find the solution G(t), which satisfies boundary 
conditions (Eil), 



G(t) = Im (g 2 e 



-z 2 t\ 



ff3 e 



-z 3 t 



(52) 



where z%, (i = 1,2,3) are the three roots of the third 
degree polynomial P(z) = (z 2 + uj 2 )(z + w c ) + jow c z. 
The constants gi are: 



.92 

93 



I ( x (w c - z 3 )(z 3 - Zl) 



Im(z 2 
(w c - z 3 ) 

\Z3~ztf 



\Z3-Zif 



(53) 



We will concentrate on the underdamped case in which 
there is only a real root (2:3) and two complex ones (z\ = 

4)- 

Given the input parameters for the problem (i.e. u , 70 
and w c ) we obtain the function G(t) and from it we eas- 
ily compute the time dependent frequency uj(t) and the 
friction coefficient "f(t). These functions are plotted in 
Figure 1 where we see how they vary on the very short 
time-scale 1 /w c (the cutoff time-scale) . From the above 
formulae we can find analytic expressions for the initial 
and final values of these functions. The time depen- 
dent frequency is initially equal to the unrenormalized 
frequency, i.e. uj 2 \q = tu 2 + 70 w c and its final value is 
equal to the renormalized frequency lu . On the other 
hand, the time dependent friction initially vanishes and 
the asymptotic value is equal to 70. 

The diffusion coefficients are temperature dependent. 
In the high temperature regime the results are well 
known: both coefficients (that start being zero) approach 
asymptotic values given by D\ = 7ofcsT and D 2 = 0. We 
studied the behavior in the zero temperature case where 
it is possible to find exact analytical expressions for both 
coefficients. The easiest way to present the results is to 
notice that the position autocorrelation function can be 
written in terms of exponential integral functions as: 

S(t) = ^Im( ff2 (e^E^i) - e-^Ei(z 2 t))) 

+ ^-g 3 (e Z3 *E 1 (z 3 i)-e- Z3 *Ei(z 3 i)) (54) 

Using this (and relating D\ and D 2 to S(t)) we obtained 
the plots of the diffusion coefficients shown in Figure 1. 



I(w) = 70W 



!(!- 



(51) 
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FIG. 1. Time dependent coefficients entering in the homo- 
geneous terms of the master equation. The environment is 
described by Drude's spectral density and the initial state 
is of zero temperature. The time dependence (in units of 
the renormalized frequency) of all coefficients show an ini- 
tial transient for times of the order of the cutoff time-scale. 
Parameters for the plot are to = 1, w c = 100, 70 = 0.01 

The time dependence of the inhomogeneous coefficients 
C\ and C2 , can also be computed in this way and the re- 
sult is shown in Figure 2. The basic feature is that both 
coefficients are exceedingly small and become neglible af- 
ter a time which is of the order of the cutoff time-scale. 
After this short initial transient, the impact of the initial 
correlations on the future evolution of the system can be 
entirely neglected. 




0.01 0.02 0.03 0.04 0.05 
time 

FIG. 2. Coefficients appearing in the inhomogeneous terms 
of the master equation. They carry the influence of the ini- 
tial correlations on the future evolution. They both become 
negligible after a time of the order of the cutoff time-scale. 
Parameters for the plot are uj = 1, w c = 100, 70 = 0.01 



V. DECOHERENCE 

Here, we will examine the time evolution of delocal- 
ized initial states analyzing the effectiveness of the de- 
coherence process. We will consider two related initial 
conditions. 



A. Superposition of two translations 



Let us first analyze the initial state: 
Po oc (e lL ° p ° + e - iL °t') P[j (e iL °P° + e ~ lL ^) (55) 



where p s is the momentum operator for the system. The 
evolution of this state is particularly simple to analyze. 
Moreover, in spite of its simplicity, this state still cap- 
tures some of the essential features present in a realistic 
"Schrodinger cat" state. We will show below that the 
conclusions we obtain for the state (55) remain qualita- 
tively correct for more realistic ones. Such state cannot 
be prepared by making a measurement on the system 
only and, roughly speaking, it represents a "superposi- 
tion of two translations": In fact, if pp is a pure coher- 
ent state for the system, p Q is a superposition of two 
coherent states each one of which is displaced by ±Lq. 
However, as p@ is a thermal equilibrium state for the 
correlated system-environment ensemble, the interpreta- 
tion of (pq) is not so transparent. The evolution of (EF 
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was first studied in |E2[ , where only the zero temperature 
case was examined using a completely different formal- 
ism. The preparation function for state ( |55| ) is the sum 
of four delta functions (arising from the matrix elements 
of displacement operators): 



A = A ++ + A__ + A+_ + A_ + 
A±± = N 2 5{x - an) S(f -n± L ) 
A ±T = N 2 5{x - Xi ± 2L Q ) 8{f - r, 



(56) 



Using these equations and the exact form of the evolu- 
tion operator we can compute the reduced density matrix 
at arbitrary times. The simple form of the preparation 
function makes most of the integrations trivial. The fi- 
nal answer can be conveniently expressed in terms of the 
Wigner function, which turns out to be: 



W = 
W±±(r,p) = 

W int (r,p) 



N 2 
2-Kqppp 

N 2 



W— + Wi, 



e 



e" A e 2q i 



(pTpq) 



cos(K r ? 



(57) 
K p p) (58) 



where the coefficients appearing in these equations are 
given by 

(59) 
(60) 
(61) 

f-p 2 - (62) 
Above, N is the normalization constant 



V'O 


= L G, po 


= f , 




- 2L 2 ° l<* 


+ Ga 6 ), 


p 2 pKp 






A 


= 2L 2 p 2 - 


K r - - 



-2L 2 p 2 



N z 



2 1 



1 



: exp( 



1 



(63) 



The interpretation of equations (|57|,^8|) is clear: The 
Wigner function is the sum of two Gaussian peaks and 
an interference term. The Gaussian peaks are centered 
around the dissipative classical trajectories determined 
by equation ( ^5|) with initial conditions r = ±Lo,p = 0. 
The spread of each Gaussian peak is constant and given 
by the equilibrium values and p 2 3 . This means that 
the individual Gaussian peaks remain "intact" along the 
evolution of the system. In this sense, they are perfect 
"pointer states" (selected by the predictability sieve cri- 
terion discussed in |27]]). Indeed, it is easy to show that 
if only one translation operator is applied in ([55]), the 
linearity of the problem implies that the entropy of the 
reduced density matrix remains constant (i.e., for a sin- 



gle Gaussian Tr{p 2 r (t)) = Tr{p r {0)) = l/^J2q 2 p 2 ). We 
remark that these Gaussian peaks are not pure states of 



the system since, due to the initial correlations, the en- 
tropy of the reduced density matrix is nonzero even at 
zero temperature. To see this it is easier to analyze the 
weak coupling limit, where 70 <C 0J o -C w c . In such case, 
the roots are Z3 1=3 w c — 70 and z-i « 7o/2 + i \/<jJq — 7q /4 
(up to terms of order 7o/u> c and u>q/w c ). Using this ex- 
pressions we can find that 



4 



— + — log — 

2 7T UJ 

1 7o , w c 
o 2 lo S — ' 



(64) 
(65) 



i.e., the equilibrium values are such that 2q 2 j i p 2 j i ^ 1. 

The interference term in ( |58| ) is centered around the 
origin (the midpoint between the two Gaussian peaks). 
The oscillatory term produces interference fringes in 
phase space (regions where the Wigner function becomes 
negative). The initial value of the coefficients is such 
that the exponential factor is unity (i.e., A\§ = and the 
fringes are oriented along the momentum direction (i.e., 
K r |o = 0, K p \o = 2Lq). When the system starts evolving 
the wavelength of the fringes becomes larger (due to the 
effect of diffusion). Therefore, the wave- vectors k p and 
K r tend to zero inducing the growth of the exponent A(t) 
and with it, the exponential suppression of interference. 
A simple expression for the exponent A(t) in terms of the 
position autocorrelation function is 



.4 



2Llp% 



S 2 



2 2 




(66) 



This formula is quite useful since it makes evident a few 
important points: First, it clearly shows that decoherence 
is produced by the same process responsible for the decay 
of the correlation function (remember that S(t) is the 
symmetric part of the position autocorrelation function) . 
Second, it also shows that the maximum attainable value 
for A(t) is 2L\p 2 p. Third, it shows that the decoherence 
time-scale is typically much shorter than any dissipative 
or dynamical time-scale in the problem. Thus, by the 
time the derivatives of the correlator S decay to half its 
initial value, the fringes are suppressed by a factor of 
order exp(— L\p 2 ^) which, for large separations, can be 
very small. We can define the decoherence time-scale 
as the time which takes for the exponent to grow to a 
number of order unity (this time-scale is clearly inversely 
proportional to Lq). An analytic expression for A(t) can 
be obtained for Drude's model both in the high and low 
temperature limits. The result (for zero temperature) is 
shown in Figure 3. From the plot we observe a sizable 
growth of the decoherence factor occurring in a rather 
short time-scale: for times of the order of the cutoff time- 
scale we have A(t = w^ 1 ) « 0.02Lq. The subsequent 
growth of decoherence is not monotonic being maximal 
when the peaks are separated in position and minimal 
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when the separation is in momentum. To estimate the 
initial behavior of A(t) it is useful to obtain an analytic 
expression for S(t) valid for short times (t <C l/ui ). In 
this case we have 



Tr(pffP) 



(70) 



S(t) « «j§ - \p 2 pt 2 + ~h{w c t), 



(67) 



where P is a projector onto a pure state of the system 
P = |* >< vj/| and the state |<3> > is itself a Schrodinger's 
cat state (i.e., a superposition of two Gaussian packets): 



h{z) 



■(e z E 1 (z) 



7o 
2ww 2 

+2(C e + \ogz)-(] 



2 K(z) 



C e )z 2 + z 2 log z ) . (68) where \&± > are such that 



where C e is Euler's constant. The function h(t) is ini- 
tially zero and grows on a short time-scale in a rather 
smooth manner. When h is small, the decoherent expo- 
nent turns out to be approximately A(t) ~ h(l + h). For 
very short times we therefore have 

7o 



A(t) « -^wir\logw e t\ + 0{w 4 c t 4 ). 



(69) 
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FIG. 3. The evolution of the exponent which is responsi- 
ble for suppressing the interference fringes. This is A(t)/L% 
for the state considered in section 5. A and A ac {t)/LQ for the 
one considered in section 5.B. Both curves are almost identi- 
cal showing a rapid decoherence for a time-scale of the order 
of w c . We also see that decoherence does not grow mono- 
tonically with time but in an oscillatory fashion, with max- 
ima when peaks have position separation and minima when 
the separation is in momentum. Parameters of the plot are 
70 = 0.01, w c = 100. 



< x|*± > = N exp 



(qTL ) 2 
2d 2 



±iP q 



(71) 



(72) 



Thus, the initial state is rather similar to the one consid- 
ered above in equation ([55]). However, state ( |70| ) could 
be prepared through (a rather idealized) measurement 
on the system. The convenience of having analyzed the 
previous example becomes evident once we notice how te- 
dious the calculations become for the initial state (|70|). In 
fact, the preparation function one can associate to state 
) fr0|) has sixteen terms (four for each of the two projection 
operators onto the Schrodinger's cat state): 



16 



\=<q i \P\q><q'\P\q! i >=VA, 



1=1 



(73) 



Each one of these sixteen terms (which we don't explic- 
itly write down) can be evolved using the exact propa- 
gator. The integrations are all Gaussian and straight- 
forward. The final result can be conveniently presented 
in terms of the Wigner function which is formed by four 
Gaussian peaks and six interference terms. Indeed, one 
has an interference term between each pair of Gaussian 
peaks (having in mind that each interference term is the 
combination of two contributions we have twelve terms 
contributing to the interference and four to the direct 
Gaussian terms; this totals sixteen accounting for all the 
terms in the preparation function). At first sight the 
existence of four Gaussian peaks may seem awkward (if 
not simply wrong) but its origin and interpretation will 
be explained below. Before, we will complete our presen- 
tation of the Wigner function. It can be finally written 
as: 

w = w a + w b + w c + w d 

+ W ab + Wbc + W cd + W bd + Wad + W ac (74) 

The Gaussian peaks are given by the following expres- 
sions: 



B. Schrodinger's cat state 

We will now consider the evolution of a delocalized 
state in a less idealized situation than the one analyzed 
above. The initial state is: 



IF, 



N (r-r-j) 2 (.P-Pi-b {r-Tj) r 

e ~^~e , (75) 



1-KO x O v 



where the index i labels the peak (i.e., i = a, 6, c, d). On 
the other hand, the interference terms Wa are 
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N _ A . _ 
e ^ e 



■no x o Tj 



*~ x e 2 "p 2 x 



x cos - bn plj )(r - r i3 ) + n plj (p - (76) 



Each interference peak is centered about the mid-point 
between the corresponding Gaussian peaks, i.e. 



n + r,- 



ft + Pi 



(77) 



Before giving any details about the many coefficients 
entering in these equations let us analyze and justify the 
existence of the four Gaussian peaks and their correspon- 
dent interferences. First, let us mention that the tra- 
jectories followed by the peaks are determined by two 
functions (tq and r\ ) in the following way: 



ro ± ri, 



(78) 



The location of the peaks in momentum is simply ob- 
tained from the velocities, i.e. pi = r^. The functions ro 
and r\ are given by: 



r (t)=L G, n(t) = 



La 



1 



s 2 
2if 



Gap, 



(79) 



The function ro corresponds to the dissipative trajec- 
tory satisfying equation (|35|). On the other hand, the 
terms proportional to r\ originate on the initial corre- 
lations (remember that for a factorizable state one has 
Q'e = 0). Moreover, we can also show that T\ vanishes 
at the initial time, when there are only two peaks in- 
stead of four and a single interference term instead of 
six. Thus, the initial correlations seem to be producing 
a rather curious effect: Each of the two Gaussian peaks 
splits in two pieces generating an "interference" term in 
between. Accordingly, the initial interference term also 
splits in four pieces. How can this be possible? To under- 
stand this we should notice that the existence of initial 
correlations implies that the evolution of each piece of 
the initial Wigner function is not independent of the ex- 
istence of the other pieces. Thus, the role of the initial 
correlations is to produce a (very short lived) force that 
kicks the center of each Gaussian away from the trajec- 
tory determined by Tq{£). However, each piece of the 
initial state produces a different kick being the net effect 
a splitting of the Gaussian peak (this can also be though 
as a type of "non-linearity" induced by the initial corre- 
lations, which enable different pieces of the initial state 
to see each other). However, for us the relevant point 
is how big is the separation between peaks and how big 
is the wavelength of the intermediate fringes. Below, we 
will show that in realistic situations (like in the Drude's 
model at zero temperature) the separation between peaks 
is much smaller than their width while the wavelength is 
always much larger than the width of the peaks. There- 
fore, in such case there are no intermediate interference 



fringes being created by the initial correlations but only 
a small distortion of the packets which are not exactly 
Gaussian (i.e., we are simply writing a deformed Gaus- 
sian as the sum of two slightly displaced Gaussian and 
an interference term). However, in other situations where 
this formalism applies (remember that the above formu- 
lae would also describe a situation where the environment 
consists of a single oscillator) this effect could be larger. 

Now let us describe the structure of the interference 
terms. It turns out that the wave-vectors associated with 
all the interference terms can also be written using only 
two functions n ro , n r \ for position and two for momen- 
tum K po and K p \. In fact, one can show that the following 
relations hold: 



where n r i and n p i are defined as 



(80) 



pli K p c d 



Thus, the only relevant functions we need to know to an- 
alyze the wavelength of the fringes are K r o . . . k p \ . The 
explicit formulae for these functions, together with the 
ones for the dispersions a x , a p and all the other parame- 
ters defining the Wigner function are listed in Appendix 
3. The expressions are more complicated than the ones 
we analyzed in the previous subsections but, again, exact 
analytic expressions can be found for Drude's model at 
zero temperature. In that case we investigated the time 
dependence of the separation between sub-peaks (i.e., 
the distance between peaks a and b, or the one between 
c and d) and we found it to be very small compared to the 
width of the peaks. Therefore no separation can be seen 
at all. This is shown in Figure 4 where we also plotted the 
ratio between the width of the peaks and the wavelength 
of the interference fringes between sub-peaks. This is 
also very small showing that no oscillations are observ- 
able. This justify our previous claim that the sub-peaks 
are only a manifestation of a small distortion in the Gaus- 
sian nature of the principal peaks. We also analyzed the 
decay of the interference terms between principal peaks 
computing the ratio A qc /Lq, which is plotted in Figure 3. 
In that Figure we see that the behavior of this quantity 
is almost identical to the decoherence factor discussed in 
the previous subsection (which indeed is much easier to 
calculate). Thus, all the conclusions regarding the effec- 
tiveness of decoherence apply equally well to the both 
subsections. 
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FIG. 4. Separation between "internal peaks" and charac- 
teristic size of the "internal fringes". This plot shows that 
the sub-peaks never separate and the interference fringes be- 
tween them are always unobservable. The sub-peaks form 
a distorted Gaussian peak which moves around as a whole. 
Parameters of the plot are: 70 = .01, Lo — 5, w c = 100. 

Finally, we also examined the time dependence of both 
the position and momentum dispersions (as well as the 
entropy of the Gaussian state, which is related to the 
product of such quantities). They are plotted in Fig- 
ure 5. Contrary to what happens with the state dis- 
cussed in the previous subsection the dispersions depend 
in time and the initially pure state gets mixed as in- 
teracts with the environment. After a sudden burst of 
entropy (which is nevertheless quite small) the entropy 
decays towards a final value which is of the order of 
the final equilibrium entropy of the subsystem (i.e., of 
the order of log(2g^p^) sa .014. It is worth mentioning 
here that the nature of the final equilibrium state can 
be examined by analyzing the evolution equation for the 
Wigner function. In the long time limit equation (|40| ) 
has a very simple form since the inhomogeneous terms 
vanish and all the time dependent coefficients approach 
asymptotic values: ^(t) — > 70, ui(t) — > u , D±(t) — > d\ 
and D2(t) — ► d%. Thus, one can show that a Gaussian 
state is the stable stationary solution provided that the 
position and momentum dispersions are = d\j 70 and 
a x = (^p + ^V^o- Thus, the role of the anomalous diffu- 
sion term is to squeeze the final equilibrium state. Its ef- 
fect at zero temperature is evident from equations ( pip5| ) 
where we see that the final state is squeezed in position 
and spread in momentum (with respect to the oscilla- 
tor's ground state). The uneven squeezing is responsible 
for the non- vanishing entropy of the equilibrium state at 
zero temperature. 



0.025 p-j 
0.02 r 

0.015 r 
0.01 r 

0.005 r 
- 



Entropy of a gaussian state 



20 40 60 80 100 

time 



0.54 - 
0.52 

0.5 M 
0.48 j- 
0.46 - 



0.54 - 
0.52 - 

0.5 
0.48 - 
0.46 - 



20 40 60 80 100 
time 



20 40 60 80 100 
time 



FIG. 5. The position and momentum spread for a sin- 
gle Gaussian peak are modified on a very short time-scale 
during which entropy is produced. Later, they both settle 
towards the equilibrium values which are = c(p — .498 



and <j„ 



Pi 



.512. Parameters of the plot are 



70 = .01,to c = 100. 



VI. SUMMARY AND CONCLUSION 

In this paper we extended previous analysis of QBM 
to a more general class of initial states containing cor- 
relations between the system and the environment. We 
derived a master equation for the reduced density ma- 
trix which is local in time but has time dependent coeffi- 
cients and inhomogeneous terms (arising from the initial 
correlations). A detailed analysis of the coefficients was 
performed for Drudc's model of an ohmic environment. 
In such case, the corrections arising from the initial cor- 
relations are very short lived decaying in the time-scale 
associated with the high frequency cutoff. However, dur- 
ing that short time they can play an interesting role. 
A point to notice is that, contrary to previous specula- 
tions H|), the diffusive coefficients of the master equa- 
tion still exhibit an initial jolt in the cutoff time-scale. 
This jolt is relevant for decoherence producing the decay 
of interference effects: at zero temperatures and for very 
small damping (70 = 0.01), the interference between two 
wave-packets separated by a distance 2Lq is suppressed 
by a factor of order exp(— 0.02 x Lq). Therefore, initial 
jolts are not associated with the absence of initial corre- 
lations. To the contrary, they are likely to be related to 
the instantaneous preparation procedure which is applied 
here. Models with non-vanishing preparation time-scale 
are currently under investigation [ p3fl . 

Finally, we would like to stress once more the simplicity 
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of the formula we obtained for the decoherence factor. In 
fact, for the delocalized initial state examined in Section 
V.A we showed that the factor suppressing interference 
fringes is exp(— A(t)j where 



A(t) = 2L 2 p 2 1 




being S(t) the position autocorrelation function (in ther- 
mal equilibrium): 



S(t) = l<{q(t),q(0)}>(3 



(82) 



This equation enables us to obtain very simple qualita- 
tive estimates on the efficiency of decoherence. In fact, it 
clearly shows that even though decoherence has the same 
physical origin as the decay of correlations, the charac- 
teristic time-scale of both processes is entirely different. 
In fact, from the above equation one can simply estimate 
the amount by which correlations must decay for deco- 
herence to occur. Thus, at the decoherence time-scale 
(when the above exponent is of order unity), the decay 
of correlations is still very small: 




(83) 



S(tdec) 



where \p = h/yj2p'j 3 is the characteristic de Broglie 

wavelength of the system in thermal equilibrium (which 
approaches h/ \J2mkBT at high temperatures and the 
spread of the ground state at low temperatures) . In Sec- 
tion V.A we used the above formula to obtain an an- 
alytic expression for the decoherence factor in Drude's 
model at zero temperature. Moreover we showed that 
even though this formula was derived under simplifying 
assumptions concerning the initial state (obtained from 
a thermal state by "superposing two translations" ) it is 
robust when applied to more realistic cases (as Figure 3 
shows). 

As a final remark, we would like to show how simply 
the usual result for decoherence time-scale in the high 
temperature limit [Q arises from equation (81): At high 
temperatures the momentum dispersion is p 2 ^ = ksT 
and the autocorrelation function decays exponentially as 
S(t) « S(0) exp(— 7 t). Thus, equation ( |8l"| ) reduces to 

cxp {-A (*)) fa cxp (-2L 2 Q p 2 (1 - cxp (-2 7o t))) 

« exp (-4 l0 tk B TL 2 ) , (84) 

which is the usual result obtained in the high temperature 
approximation. We believe our equation will be useful for 
estimating the time-scale of decoherence in many other 
systems where the behavior of position autocorrelation 
function is well known. 
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APPENDIX A: COEFFICIENTS DETERMINING 
THE EVOLUTION OPERATOR 

The following kernels appear in the Generalized Influ- 
ence Functional: The ordinary noise kernel (extended to 
the complex plane) v{z) and the kernel Ki{z) are: 



- +oo 

(s - it) = - ^2 g n (s) exp(ii/„r) 



n— — oo 



k 2 (s - it) = i ^ f n {s) exp(ii/„r) 



TL— — OO 



k(s — IT) — v(s — IT) + tK2(s — IT) 



(Al) 



(A2) 
(A3) 



where the functions g n (s) and f n (s) are defined in terms 
of the spectral density as {y n — 2irn/ (3 are the Matsubara 
frequencies) : 



9n(s)= I —I(u>) 2 2l f 2 cos(^s), (A4) 

7T UJ Z + V± 



f n (s)= [ ^/H^^sinM. (A5) 

J 7T iO A + V* 



The Euclidean integral brings another contribution to the 
noise which turns out to be 



,!W I{u)-^ 



(A6) 
(A7) 



The explicit form of the coefficients that define the 
evolution operator is: 

a = [2 7 rG(t)(2 7 r ? |)i]- 1 (A8) 

a ^ = w> a ^ v i (A9) 

G 1 G 2 

"3 = ^, «4 = -q. a 5 = -[—-G] (A10) 



a e = ± / ds Ci(s) vi(s) 

7 J 2 

o 



(All) 
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as = ± / ds C i (s) V2 (s) 



(A12) 



t t 



1 



a 10 = — j ds J du R(s,u) vi (s) vi (s) (A13) 

o o 
t t 

an = J ds J du R(s,u) vi(s) V2(s) (A14) 



where the constants g 2 , and p| are expressed (in terms of 
the spectral density and the initial temperature) as 



%1 



n— — oo 



n— — oo 

2 i ..2 i \-l 



U n = {UJ Q +V n +XnY 



(A15) 

(A16) 
(A17) 



The auxiliary functions v±(s), V2(s), Ci(s) and C*2(s) ap- 
pearing in the above expressions are: 



wi(s) 



Gft - s) 

cm < 



«a(s) = G(i-a) 



W) 



+ OO 



n=— oo 
j +oo 

C 2 (s) = - u„ ^„ /„(s) 



while the effective noise kernel R(s, u) is 

R(s , u) — v(s — u) + V C orr{s, u) 

v corr {s,u) = -q 2 C 1 (s)C 1 {s) + 



(A18) 



£ E U n[9n{s)g n (u) ~ fn(s)fn(u)}. (A19) 



APPENDIX B: AUTOCORRELATION 
FUNCTION 



-T 

OP, ^ 



t-t' 



ds coshz/„s 



2/3 — jo 
x {G(t - t' - s) - G(t' - t + s)) 

where G is the Laplace transform of G(t). 



(Bl) 



APPENDIX C: THE WIGNER FUNCTION 

The following formulae determine the temporal depen- 
dence of the coefficients determining the Wigner function 
for a Schrodinger cat initial state: 

G 2 



2o x 2 = 6 2 G 2 
+ G 2 



5 2 

4a io + 



2q 



Ox 



G 2 



G' 



2 

<T X K r l 



4«?2 + 1 f 1 

^ + GS 2 P 
PoGae 



2 " r 4<5 2 



L Gai 



2<l 



Ox d 

cr„Kp° = — cr r K™ 



A. 



2p 2 <5 2 



(CI) 
(C2) 

(C3) 

(C4) 
(C5) 

(C6) 



2^ 2 + l £ + 1 



£ 2 



j8 



2 ^.2 
— ^— (ft po ± Kpi) 2 — (K ro ± K r i) 2 (C7) 



.4, 



1 



* + 1 2p 2 <P + i 
tip 

-r 2„ 2 _ 2,, 2 

Op fupi /VI 



r 2 



(C8) 

T 2 a 2 k 2 rr 2 k 2 

A — J_ P 2 A 2 _ p po x ro (C<Q\ 

A n - -p + ^° 6 2 2 1 j 



The following formula (which is proved in jllj] using 
Laplace transform techniques) enable us to obtain simple 
relations between the auxiliary function U(t,t') and the 
position autocorrelation function S(t): 

U(t, t') = P lG(t)G(t') + S(t)G(t') + S(t')G{£) 

+ S{t)S y ) G(\v n \)coshv n (t - t 1 ) 



4 
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